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ABSTRACT 

Aims. The Luminous Blue Variable (LBV) R71 is currently undergoing an eruption, which differs photometrically and spectroscopi- 
cally from its last outburst in the 1970s. Valuable information on the physics of LBV eruptions can be gained by analyzing the spectral 
evolution during this eruption and by comparing R71's present appearance to its previous outburst and its quiescent state. 
Methods. An ongoing monitoring program with VLT/X-shooter will secure key spectral data ranging from visual to near-infrared 
wavelengths. Here we present the first spectra obtained in 2012 and compare them to archival VLT/UVES and MPG/ESO- 
2.2m/FEROS spectra from 2002 to 2011. The discussed data include pre-eruption spectra in 2002 and 2005, a spectrum of the 
transitionary phase between quiescent and eruptive state in 2007, and spectra of the eruption in 2011-2012. Information on R71's 
1970s outburst is taken from the literature. 

Results. The 2011-2012 spectra are dominated by strong neutral and singly ionized metal absorption lines likely formed in a large 
"pseudo-photosphere." We find an unusually low apparent temperature of R71 of only r e ff 2oi2 ~ 6 650 K; the star resembles a late F 
supergiant. R71's visual lightcurve had a maximum in 2012 with mv,20i2 ~ 8.7 mag. Given the uncertainty in the extinction towards 
R71, this corresponds to a bolometric luminosity of Mb o i,20l2 ~ _ 9.8 mag to -10.3 mag. R71's 2011-2012 spectra do not show 
HI and Fell P Cyg profiles, which were present during its last outburst in the 1970s and which are normally observed during LBV 
outbursts. Low-excitation forbidden emission lines and Fe I P Cyg-like profiles from a slowly expanding nebula became apparent in 
late 2012. These lines originate likely in the rarefied region above the pseudo-photosphere up to 13 AU from the star. 
Conclusions. The rise in R71's visual magnitude and the low apparent temperature of its pseudo-photosphere during the current 
eruption are unprecedented for this star. R71 most likely increased its bolometric luminosity by AMy = 0.4-1.3 mag compared 
to its quiescent state. The very low temperature of its pseudo-photosphere implies a very high-mass loss rate on the order of 
^R7 1,20 12 ~ 5xl0~ 4 M Q yr~' compared to M qu j escence ~ 3xlO~ 7 M G yr -1 . The apparent radius increased by a factor of 5 to about 
500 R Q . No fast-moving material indicative of an explosion is observed. The changes in R71's photometry and spectrum are thus 
likely consequences of a tremendously increased wind density, which led to the formation of a pseudo-photosphere. 

Key words. Stars: variables: S Doradus - Stars: winds, outflows - Stars: mass-loss - Stars: individual: R71 



• i— i ' 1 . Introduction LBVs, also known as S Doradus variables, are evolved mas- 

/\ , sive stars that exhi bit a particu l ar type of instability which is 

^ . The LBV R7 1 (= HDE 269006) in the Large Magellanic Cloud not yet understood dContill 1984 fl997l |£phrevs & Davidson 

_ C5 _■ (LMC) is currently undergoing an eruption that started in 2005 m4 ^ No ta & Lamerslll997l and references therein). They rep- 

dGamen et al. | 2009|) . In 2012 the star had reached unprecedented resent a brief but critical phase in massive star evo i ution be - 

visual brightness accompanied by remarkable variations in its caufje several solar masses can be expelled during this stage . 

optical spectrum ( |Gamen et al. || 2012|) . R71's current eruption is Most of the f undamen tal questions about the physical cause 

characterized by more extreme parameters than its outburst in of their instability are sti n unsolved and they may fall into 

the 1970s and may provide us with important insights into the different cate gories, i.e., classical LBV outbursts and giant 

LBV phenomenon. eruptions. Outbursts with visual magnitude variations of 1- 

2 mag and constant bolometric luminosity are commonly re- 

Send offprint requests to: A. Mehner, e-mail: amehner@eso.org f err ed to as classical LBV outbursts. During giant eruptions 

* Based on observations collected at ESO's Very Large Telescope un- the visual magnitude increases by more than 2 mag and 

der Prog-IDs: 69.D-0390(D) and 289.D-5040(A) and at the MPG/ESO the bolometric luminosity likely increases. Historical exam- 

2.2-m Telescope under Prog-IDs: 076.D-0609(A), 078.D-0790(B), p i es of g iant er uptions are P Cygn i in the 17th century (e.g., 

086.D-0997(A), and 087.D-0946(A). de GroQ f 19g8 Lamerfj & de ^ 1992)> Cal in th / ^ 

This work was co-funded under the Marie Curie Actions of the / T It~, j o tt u — l N nml Unnl a 

„ . . (see the reviews by Davidson & Humphreys 1997; 2012 and 

European Commission (FP7-COFUND). J r 
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references therein), V12 in NGC 2403 in the 1950s (e.g., 
iTammann & Sandagdll968l ISmith et alJl200ll). and SN1961V 
in NG C 1 058 (e.g..lGoodrich et al.lll989l IVan Dvk & Mathesonl 
l2012t cf. [Kochane k et alJl201 lh . The most promising explana- 
tion for the LBV instability mechanism involve radiation pres- 
sure instabilities, but also turbulent pressure instabilities, vibra- 
tions and dynamical instabilities, and bi narity cannot be entirely 
ruled out (Hum phreys & Davidsonl 19941 and references therein). 
It is unknown what the relative roles of classical LBV outbursts 
and giant eruptions are and which of the two events are more 
important in terms of total mass lost. It is also not established if 
giant eruptions result from unusual circumstances and due to a 
different instability mechanism or if they are only extreme cases 
of classical LBV outbursts. 

LBVs may be critical to our understanding of massive 
star evolution. Their instability could be responsible for the 
empirically found upper luminosity boundary (also known as 
Humphreys-Davidson li mit) above and to the right of whic h 
no supergiants are found (Humphreys & Davidso nlll979Hl984h . 
The commonly accepted idea is that very high-mass loss pre- 
vents stars above 40-50 M from becoming red supergiants 
(RSGs) but instead drives them towards higher temperatures. 
LBVs are located in the region of concern in the HR-diagram and 
their instability leads to high-mass loss episodes. The realization 
that winds of O and Wolf-Rayet (WR ) stars are clumpy/porous 
( Full erton et al. 2006, Pu is et al.l200 8) has reduced their inferred 
mass loss rates by a factor of ~10. This further increased the in- 
terest in the LBV phase, since high-mass loss episodes may be 
required for massive st ars to evolve directly to the WR stage 
(ISmith & Owockill2006h . 

LBVs are g enerally considered to be stars in transition to 



the WR stage dMaederl Il983l iHumphrevs & Davidsonl 1 1994 
Lang er et al.l ll994). but observations of LBV-like eruptions im- 
mediately prior to supernova (SN) explosions and also recent 
stellar evolution modeling may lead to a revision of this inter- 
pretation. For several exceedingly luminous SNe it has been 
suggested that shock interaction with dense and opaque cir- 
cumstellar material cou ld account for the excess luminosity 
(ISmith et alJ 12007: 2008). The progenitor stars likely expelled 
extensive circumstellar ejecta durin g LBV outbursts or giant 
eruptions. lOuataert & Shiodel (120 1 2h proposed an alternative ex- 
planation for large stellar mass-loss rates just prior to core col- 
lapse. During core Ne and O fusion in massive stars the fusion 
luminosity exceeds the star's Eddington luminosity, which drives 
convection and in turn can excite internal gravity waves. The 
gravity waves convert into sound waves as they propagate to- 
wards the stellar surface and the subsequent dissipation of the 
sound waves can unbind several solar masses of the stellar enve- 
lope. 

In a few cases LBV-like e ruptions have been observed prior 
toaS N explosion; SN2006ic jFolev et al 1120071 iPastorello et alJ 
20071) . SN2005gl dGal-Yam et alJ 120071 jGal-Yam & Leonard! 
2009), and SN2 009ip dMauerhan et al.l 120121) but compare 



Dwarkadas (201 1) who argued that only in the case of SN2005gl 



there is substantial evid ence for an LBV p rogenitor. In the re- 
cent case of SN2009ip. fSmith et al.l (|2010) identified the pro- 
genitor as a likely LBV based on its photometry but its mul- 
tiple pre-SN outbursts were most likely due to pulsational pair 
instability, outbursts associated with late stages of nuclear burn- 
ing, or close periastron passages of a companion star rather than 
wind-driven LBV outbursts because their time sc ales were much 
shorter than for mos t other LBV outbursts dMauerhan et al.l 
[2012b . ISoker & Kashil d20H proposed that SN2009ip's pre-SN 
outbursts were due to binary interaction and that its final ex- 



plosion was a mergerburst and not a true SN. lOfek et al.l (120 131) 
claimed that a mass loss event of the progenitor of SN2010mc 
is consistent with a wave-driven pulsation model and is causally 
connected to its SN explosion only 40 days after. 

The link between SNe and L BVs had no t been supported by 
stellar evolutionary models until iGroh et all (1201 3l) reported on 
their theoretical finding that single rotating stars with initial mass 
in the range of 20-25 M have spectra similar to LBVs before 
exploding as SNe. This result sets the theoretical ground for low- 
luminosity LBVs, such as possibly R71 (but see Section @), to 
be the endpoints of stellar evolution and may also strength en the 
identification of the progenitor of SN1987A as an LBV (iSmithl 
12001 . 

Important clues towards an understanding of the LBV phe- 
nomenon can be gained from their variabilities that occur on dif- 
ferent magnitude and time scales (IHump hrevs & Davidsonl 19941 
and references therein). R71 is excellent for a case study. It has 
been well-observed over the last decades and its location in the 
LMC is fortunate because its distance and reddening are well 
constrained. During its quiescent phase R71 has a hot supergiant 
spectrum with strong H I, He I, Fe II, and [Fe II] lines and P Cyg 
profiles. It was cla ssified as a B2.5 Ieq star with my ~ 10.9 mag 
(Fea st et all [T960). It has a typical su pergiant mass loss rate 
of ~3xl0~ 7 M Q yr" 1 dWoif et al] fl98lT> . The star shows mi- 
crovariations of Amy ~ . 1 mag on time scales of 14-100 days 
dvan Genderen et al.ll985t[l988[[T997llLamers et al.ll998l) . R71 
had a classical L BV outburst in 1970-1977 when it reached 
my ~ 9.9 mag dThackeravl 11974 IWoiflll975l Ivan Gendereiil 
U979llWolf et al.ll981l.lvan Genderen! 19821) . During the outburst 
its mass l oss rate increased b y a factor of about 100 t o ~5xl0~ 5 
Mq yr" 1 dWolf et al.lll98ll cf. lVink & de Koteril2002l who found 
that R71's mass loss rate increased only by a factor of 3-5). 

R71's current eruption started in 2005. In 2012 the star 
had reached an unprecedented visual magnitude of my ~ 
8.7 mag and is currently the visually brightest star in the 
LMC ( |http://www aavso.org} cf. lFeast etai](l960l) . At the end 
of 2012 a decline in R71's visual lightcurve commenced while 
the R-band l ightcu r ve sta yed bright (http://www.aavso.org). 
iGamen et ail d2009h 1201 2|) found that the spectrum in 2008 
February resembled an early-A supergiant, in 2009 August an 
extreme early-F hypergiant, and in 2012 that of an early-G su- 
pergiant. 

In this paper we discuss R7 1 's current eruption in relation 
to its previous outburst in the 1970s and to its quiescent state. 
Because of the rarity of LBV outbursts and their potential impor- 
tance in massive star evolution this large-amplitude event calls 
for special attention. In Section|2]we describe the observations. 
In Section[3]we analyze the changes of selected spectral features 
in detail. In Section|4]we discuss the results and consider R71's 
classification as an underluminous or classical LBV. In Section 
[5] we summarize our conclusions. 



2. Observations and Data Analysis 

In 2012 August, September, and November we obtained the 
first X-shooter spectra of R71 of a long-term monitoring pro- 
gram at the Very Large Telescope (VLT). X-shooter is a 
medium-resolution echelle spectrograph that simultaneously ob- 
serves with 3 arms covering the wavelength region from 3 000- 
24 800 A dVernet et al.l 1201 ll) . Spectra were obtained with the 
narrowest available slits of 075 in the UVB arm, 074 in the VIS 
arm, and 074 in the NIR arm yielding spectral resolving powers 
of R ~ 10000-18 000. In addition, spectra were obtained with 
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Table 1. Journal of observations. 



Date 


MJD 


Instrument 


Wavelength Range 


Total Exposure Time 




(days) 




(A) 


(s) 


2002-07-23 


52478.4 


UVES 


3 260-4450 


360 


2002-07-23 


52478.4 


UVES 


4 580-6 680 


3x110 


2002-08-03 


52489.4 


UVES 


3 260-4 450 


360 


2002-08-03 


52489.4 


UVES 


4 580-6 680 


3x110 


2002-08-12 


52498.4 


UVES 


3 260-4 450 


360 


2002-08-12 


52498.4 


UVES 


4 580-6 680 


3x110 


2002-09-04 


52521.4 


UVES 


3 260-4450 


360 


2002-09-04 


52521.4 


UVES 


4 580-6 680 


3x110 


2002-09-10 


52527.4 


UVES 


3 260-4450 


360 


2002-09-10 


52527.4 


UVES 


4 580-6 680 


3x110 


2002-09-18 


52535.3 


UVES 


3 260-4 450 


360 


2002-09-18 


52535.3 


UVES 


4 580-6 680 


3x110 


2002-09-22 


52539.3 


UVES 


3 260-4 450 


360 


2002-09-22 


52539.3 


UVES 


4 580-6680 


3x110 


2002-09-30 


52547.4 


UVES 


3 260-4 450 


360 


2002-09-30 


52547.4 


UVES 


4 580-6 680 


3x110 


2002-10-06 


52553.4 


UVES 


3 260-4 450 


360 


2002-10-06 


52553.4 


UVES 


4 580-6 680 


3x110 


2003-12-12 


con i £ i 

53716.2 


I I 1 W li 

rEKOS 


n zr\r\ r\ inn 

3 500-9 200 


2x450 


200/ -02-22 


Z A 1 CO O 


rbKU>S 


i zr\r\ c\ inn 
3 500—9 200 


3x500 


ZU1 1-lO-ZU 


JJ041 .U 




d jUU— y ZUU 


loUU 


ZU1 1-U J- LJ 


JJ07 / .U 




3 jUU— y ZUU 


loUU 


ZUlZ-Uo-ZO 


JOlOJ.3 


X-shooter 


1 n iac\ ia cnn 
1U Z4U— Z4 oUU 


jXOU , ZX1U 


2012-08-26 


56165 3 


X-shooter 


3 000-5 595 


3x120" , 2x20* 


2012-08-26 


56165.3 


X-shooter 


5 595-10 240 


3x120"! 2x20* 


2012-09-04 


56174.3 


X-shooter 


10 240-24 800 


3x60" ,2x10* 


2012-09-04 


56174.3 


X-shooter 


3 000-5 595 


3x120" , 2x20* 


2012-09-04 


56174.3 


X-shooter 


5 595-10 240 


3x120" , 2x20* 


2012-11-01 


56232.2 


X-shooter 


10 240-24 800 


3x60" ,2x10* 


2012-11-01 


56232.2 


X-shooter 


3 000-5 595 


3x120" , 2x20* 


2012-11-01 


56232.2 


X-shooter 


5 595-10 240 


3x120" , 2x20* 



Notes. (fl) Slit width: 07 5/0".' 4/07 4 for the UVB/VIS/NIR arms. <*> Slit width: 5"/5"/5" for the UVB/VIS/NIR arms. 




4250 



4500 



4750 



4800 



4850 



4900 



4950 



5000 



Wavelength (A) 



Fig. 1. Spectral variations between R71's quiescent state in 2002 July and its eruptive state in 2012 November. The 2002 UVES 
spectrum shows strong [Fell] emission and Fell P Cyg profiles (e.g., at AA92A A). H/3 has a P Cyg profile, while higher Balmer 
lines appear only in absorption. The 2012 X-shooter spectrum is dominated by strong metal absorption lines. (The continuum is 
normalized to unity. The wavelength region /M4500^-750 A is not displayed because of a gap in wavelength between the UVES 
blue and red settings.) 
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the 5" slits in all 3 arms to achieve an absolute flux calibration 
but also to investigate the spectral energy distribution. 

We retrieved archival spectra of R71 obtained with 
VLT/UVES in 2002 and with MPG/ESO-2.2m/FEROS in 2005, 
2007, and 2011. The UVES spectra in 2002 were obtained 
during R71's quiescent state. They cover the wavelength re- 
gion from 3 000-6 800 A with spectral resolving power of R ~ 
40 000. The FEROS spectra cover the wavelength region from 
3 500-9 200 A with spectral resolving power of R ~ 48 000 and 
were obtained throughout the current outburst starting at its on- 
set in 2005. 

A journal of observations is listed in Table [TJ Each data set 
was reduced with the corresponding ESO pipeline (X-shooter 
pipeline version 1.5.0, UVES pipeline version 5.0.17, FEROS 
pipeline version 1 .60). The spectra were corrected for the instru- 
ment response and atmospheric extinction but not for interstellar 
extinction and telluric lines. 

We used the 2002 UVES spectra to determine R71's sys- 
temic velocity. During the quiescent state the permitted lines 
show P Cyg profiles and we therefore used the [Fell] emis- 
sion lines at ^4287,4414,4416,4452,4458,4489,4728,4890 A, 
which originate in the outer stellar wind. We found 
v sys = 192+3 km s _1 , in good agreement with the literature val- 
ues v svs = 195 km s _1 dThackeravll 1 9741) and v sys = 193 km s _I 
(IWolfetalJ [l98lh . which were also determined from [Fell] 
lines. All radial velocities stated below are relative to the sys- 
temic velocity. 

3. Results 

The 2012 X-shooter sp ectra confirm the dramatic spec- 
tral changes reported in iGamen etaL I d20L2l) . However, the 
peak visual magnitude appears to have been lower (see 
http://www.aavso.org) as stated by Gamen et al. and we do not 
confirm that R71 currently resembles an early G-supergiant. 
Figure[TJcompares a UVES spectrum obtained in 2002 July with 
an X-shooter spectrum obtained in 2012 November. Most emis- 
sion lines have disappeared and the spectrum is dominated by 
strong neutral and singly ionized metal absorption lines (e.g., 
Fe I, Fe II, Si I, Si II, S I, N I, and O I) indicative of a cold dense 
pseudo-photosphere, which develops during an LBV o utburst 
( Leitherer et al 1 119851 lAt)Denzellerlll986l lDavidsonlll987l) . 

We analyzed several spectral features between 2002-2012. 
The UVES spectra in 2002 were obtained during the maxi- 
mum of a microvariation (my ~ 10.7 ma g compared to my ~ 
10.9 mag in quiescence, see Figure 1 in lSzczvgiel et al.ll2010l) . 
Small variations in the line profiles occur during 2002 but the 
spectra nevertheless closely resemble R71's appearance in qui- 
escent phaseQ At the time the 2005 FEROS spectrum was ob- 
tained, the current eruption had already commenced (my ~ 
10.7 mag). The 2002 July and 2005 December spectra show no 
major differences. The 2007 FEROS spectrum appears to be a 
snapshot of the transitionary phase between the quiescent and 
the eruptive state with my ~ 10.3 mag. The 2011 FEROS and 
2012 X-shooter spectra show R71 in eruption when the visual 
magnitude reached up to my ~ 8.7 mag (http://www.aavso.org 
cf. lGamen etai]l2012l who observed a maximum brightness of 
my ~ 8.3 mag). 

1 No significant line variations hinting at a potential companion is 
observed during the 3 months of UVES coverage. Comparison of the 
UVES spectra to the FEROS and X-shooter spectra with respect to line 
variations due to binarity is inhibited because of the spectral changes 
inflicted by the eruption. 



Figure|2]shows a time series of Ha and H/J from 2002-2012. 
Ha changed from a P Cyg profile with broad emission wings 
extending up to v em ~ ±850 km s _1 in 2002-2005 to a promi- 
nent double-peaked symmetric profile in 2011/2012. The max- 
imum depth of the P Cyg absorption in 2002 July is at v a b s = 
-118+12 km s _1 and the blue edge is at v e d ge = -186+5 km s . 
We use the blue edge of the Ha P Cyg absorption as a proxy for 
the terminal velocity but since the absorption is not well-defined 
and filled in with extra emission this measurement gives only a 
lower limit. The profile shows two additional absorption features 
at v a b Si i = -37±2 km s _1 and at v a b s ,2 = 2±6 km s . The for- 
mer may be due to an expanding shell, while the latter originates 
most likely from the photosphere. In spectra obtained between 
2002 August and October we find v e d ge = -131+20 km s _1 . In 
2007 extra emission compared to the continuum is observed. The 
P Cyg absorption has disappeared as has the absorption compo- 
nent at v a bs,i 37 km s _1 . The weak absorption component at 

system velocity strengthened, most likely due to the formation 
of a cooler pseudo-photosphere. In 2011/2012 this absorption 
component dominates the Ha profile. It has a full width at half 
maximum (FWHM) of 44± 1 km s _1 . There was no indication of 
this absorption component during the previous outburst but this 
absorption feature was first observed in the 1990s. Weak broad 
Ha emission is still present and a weak absorption component 
can be observed with a blue edge at v e d ge = -199+9 km s . 

We found comparable values for the blue edge of the Ha ab- 
sorption during quiescence and eruption. However, because the 
Ha absorption depends on the density and ionization structure 
of the wind, and the star is currently hidden beneath a pseudo- 
photosphere, we cannot state for certain whether the wind ve- 
locity stayed unchanged during the current outburst. Our val- 
ues of v e dg e 131 km s _1 to -199 km s _1 are comparable to 

the va l ue v edre = -158 k m s 1 in 1984 found bv IStahl & Woil 
(1986). W olf et all (11981) found a maximum outflow velocity of 
-127 km s" 1 in UV Fe II lines. 

H/3 behaves similarly to Ha but does not show prominent 
line emission in 2011/2012. Higher Balmer lines are always in 
absorption. Paschen lines and Bry, which were likely in emission 
during quiescent state, are in absorption in 20120 The absence 
of hydrogen emission (apart from Ha emission) indicates that 
R7 1 's pseudo-photosphere is much cooler and denser during the 
current eruption than during i ts 1970s outburst when Ha to H6 
showed strong P Cyg profiles (IWolf et al.lll981l) . 

Figure [3] compares Call lines in 2005 and in 2012. These 
lines provide information on the (changing) conditions of the 
pseudo-photosphere and the circumstellar material. Call has a 
low ionization potential and the broad and very deep Ca II H and 
K absorption lines observed in 2011 and 2012 are therefore an 
indicator of the low temperatures in R71's pseudo-photosphere. 
The forbidden [Call] AA729 1,7324 lines are in emission and 
the near-infrared triplet Call /U8498,8542,8662 is in absorp- 
tion. The Call near-infrared triplet lines are usually observed 
in emission in warm hypergiants. They are formed by radiative 
de-excitation from the upper level of the Ca II H and K absorp- 
tion, which leaves the Ca + ions in the upper level for the [Ca II] 
/U7291.7324 linesQ We investigated the line variations in our 



2 No archival near-infrared spectra of R7 1 are available and since this 
paper focuses on the comparison between R71's current eruption with 
its quiescent phase and its previous outburst we will only mention a few 
near-infrared features. R71's near-infrared spectrum during the current 
eruption will be discussed in more detail in a forthcoming paper. 

3 Alternatively. iPrieto et all ( 120091) and lSmith etall d2009l) have sug- 
gested that the [Call] AA729 1,7324 emission lines may be due to dust 
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-1000 -500 500 1000 -500 -250 250 500 
Velocity (km s" 1 ) 

Fig. 2. Ha and HyS in R71 from 2002-2012. The Ha P Cyg pro- 
file changed to an inverse double-peaked symmetric profile dur- 
ing the eruptive state. The absorption component close to system 
velocity now dominates the profile. H/3 behaves similarly to Ha. 
In the latest spectrum from 2012 November, the Fe I /16546 emis- 
sion line becomes apparent towards the blue of Ha. The strong 
absorption features to the red and blue of H/3 in 201 1-2012 are 
blends of Fel, Nil, Til, Till absorption lines. (The continuum 
is normalized to unity, velocities are in R71's restframe.) 



2012 spectra in order to explain the puzzling combination of 
Call near-infrared triplet absorption and [Call] AA729 1,7324 
emission observed in R7 1 . We find that the Ca II near-infrared 
triplet lines show the same residuals as the Ha emission, while 
pure photospheric lines vary little. The Ca II triplet is likely dom- 
inated by absorption but filled in by an emission component from 
the circumstellar material that may be sufficient to populate the 
upper level of the [CaE] ^729 1,7324 lines. 

Other spectral lines confirm the very low temperature of 
R71's pseudo-photosphere. In the previous outburst almost all 
lines (in particular the Balmer lines, Fell, and Till) showed 
classical P Cyg profiles with a deep absorption component 
(IWolf et alJll98lh . which is not the case for the current erup- 
tion. [Fell] emission lines have disappeared. Mgll /14481 
absorption, indicative of low temperatures, is strong. He I 
/L14388,447 1,5876,6678 absorption at system velocity in pre- 
emption spectra have disappeared in 2011/2012. We have no 
comparison spectra for the near-infrared He I /110830 line dur- 
ing the quiescent state and thus cannot conclusively identify the 
absorption feature and the small emission observed in 2012 in 
the concerning wavelength region as He I. More likely no He I 
absorption and emission lines are present at the current time. 
During the 1970s outburst the spectrum showed He I absorption 
(Thack eravl[T974l) . again indicating lower temperatures during 
the current eruption. The Nal /L15890,5896 lines show three 



destruction. However, it is unclear what would cause the destruction of 
dust during this eruption. 




| 7280 7290 7300 7310 7320 7330 7340 



(c) Call 8498 Call 8542 Ca II 8662 




8450 8500 8550 8600 8650 
Wavelength (A) 



Fig. 3. Call lines in 2005 and 2012; (a) the Call H and K 
lines show very broad and deep absorption, (b) the [Call] 
zl/17291,7324 lines are in emission, and (c) the Call near- 
infrared lines at /L18498, 8542, 8662 A are in absorption. The 
Paschen lines /1/18502, 8545, 8665 A can be seen in absorption 
just to the right of the Call triplet lines. (The continuum is nor- 
malized to unity.) 



absorption components: 1.) the Galactic interstellar absorption 
with v = -177+1 km s _1 relative to R71, 2.) an absorption com- 
ponent with v = 13+1 km s , which likely is the interstellar 
LMC feature, and 3.) a component with v = -29+1 km s , 
which is conceivably formed in the cool pseudo-photosphere 
and the circumstellar material based on its line strength varia- 
tions. The velocity of this component stayed constant compared 
to pre-outburst, but the line became stronger. 

Between 2012 August and November Fe I emission lines ap- 
peared redwards of their absorption components (producing a P 
Cyg-like profile). In Figure [2] e.g., the Fel /16546 line can be 
seen in the l atest spectrum (at abo ut -750 km s _1 from Ha). 
According to lKovtvukh et all d201 ll) the observed Fe I A643Q to 
Fel /16421 ratio of ~2, i.e., different from about unity, would 
indicate a non-thermal excitation mechanism. Other emission 
lines characteristic of a very low-excitation "nebula" (excitation 
energy on the order of 4 eV and less) appeared and became 
stronger. We tentatively identify narrow emission from, e.g., 
[OI] ^6300,6364,6392, [Nil] A7394, and [SI] AA7725,10821. 
These lines indicate the presence of a neutral nebula, likely 
ejected during the current eruption and possibly the rarefied 
region above the pseudo-photosphere. The Fe I emission lines, 
which appeared in late 2012, may originate from the same re- 
gion. Classical nebular lines from an ionized nebula such as 
[Nil] /U6548, 6583, which were observed in R71's spectrum in 
1984 (IStahl & Wolf! 1 19861) . are not present. The expansion ve- 
locity of the nebula is on the order of a few km s _I and thus 
small compared to the stellar wind velocity. The eruption started 
about 6 years ago and, if we assume an expansion velocity of 
Vemption ~ 10 km s _1 , the ejected material is now reaching a dis- 
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tance of about 1 3 AU or 3 000 R e from the star. At a distance of 
the LMC, this corresponds to about 0.3 mas. 

Photometrically and spectroscopically, R71's current erup- 
tion is very different compared to its previous outburst. Its vi- 
sual lightcurve reached a peak magnitude in 2012 that was about 
1 mag brighter than during the 1970s outburst and R71's current 
spectrum implies a much cooler pseudo-photosphere. 

4. Discussion 

R71 was well-observed over the last decades covering its qui- 
escent phases and its 1970-1977 classical LBV outburst. This 
gives us a unique opportunity to compare the current eruption 
to its previous behavior. Table |2]gives an overview of R71's pa- 
rameters found for its quiescent state and its 1970s outburst by 
different authors, and parameters estimated in this paper for its 
current eruption. 

There is ambiguity in the lit erature regarding th e classifica- 
tion of R71 as a classical LBV (Lennon etaTII 1993b or as a less 
luminous LBV (IWolf et all 1 1981 7 Reasons are the uncertainty 
in the extinction towards R71 and the bolometric correction, see 
Table [2] To interpret R7 1 's behavior it is important to know to 
which of the two groups it belongs because they diff er in mass 
and evolutionary path dHumphrevs & Davidsonll9 94). Classical 
LBVs have M\, Q \ < -9.6 mag and have very likely not been 
RSGs. Less luminous LBVs have Mboi = -8 mag to -9 mag, 
lower temperatures, smaller amplitudes of their outbursts, and 
lower mass loss rates. They have probably been RSGs. 

Support for the classification of R7 1 as a less luminous LBV 
comes from the fact that its visual light and temperature vari- 
ation during the 1970s outburst showed much sma ller varia- 
tions t han is observed f or the more luminous LBVs (iThack erav 
11974 IWolf et alJ 119811) . Further evidence in favor of a RSG 
phase comes from infrared studies. A dust shell around R71 
may have been produced i n a short RSG ph ase dVoors et al.l 
1999, Ivan Loon et all 1201 Ol lBoveretal.ll2010l). If the period- 
luminosity relation for LBVs found bv lStothers & Chinl |l995) 
holds for R71, then the timescale of about 40 years between the 
current eruption and its 1970s outburst implies M\, Q \ ~ -9 mag 
(only 20 years between outb ursts would have be en expected in 
case of Mboi ~ -10 mas). ILennon et al .] (119931) . however, ar- 
gued that the extinction towards R71 is anomalous and greater 
than previously thought and thus found a much higher luminos- 
ity and temperature, more consistent with a spectral type of B2.5 
Ieq. R71 enters the region of classical LB Vs in the HR-d i agram , 
if we adopt the extinction value found bv lLennon et alJ 0993). 
Its 1970s outburst may have been unusually weak compared to 
other classical LBV outbursts. 

R71 does not appear to lie in a region with anom alous 
extinction, see, e.g., Figure 1 1 in llmara & Blitzl (120071) and 
Figure 5 in Dohashi et al. (20081). The empirical relation between 
sodium absorption and dust extinction found bv lPoznanski et all 
(2012) and Ry — 3.1 implies an interstellar extinction of 
Ay = 0.09+0.03 mag t owards R71. This value is likely too 
small. iBover et al.l(l2010|). e.g., assumed Ay = 0.4 mag based on 



the extinction map by Schle geTet alJ (Q998). A potential larger 
extinction towards R71 would likely be of circumstellar nature 
and may be caused by dust formed in carbon deficient ejecta 
(ILennon et alJll993l) . An upper limit to the extinction towards 
R71 of Av,20i2 < 0.7+0.4 mag can be inferred from the ab- 
sence of the A 10780 absorption line for which a relati on between 
the lin e equivalent width and E(B-V) was found bv iGroh et al.l 
(2007). While this upper limit does not discriminate against any 
of the extinction values found in the literature, it implies that 



the current extinction value is certainly not much higher than the 
pre-emption value. The use of other diffuse interstellar medium 
lines and their correlation with extinction is inhibited because 
line identification is made impossible by many small emission 
and absorption lines even in R71's quiescent state. R71 has likely 
an extinction value of 0. 1 < Ay < 1.1 mag and thus lies close to 
the boundary between the less luminous and the classical LBVs 
during its quiescent state and a previous short RSG phase cannot 
entirely be ruled out. 

The observed spectral changes during R7 1 's current eruption 
are expected in an outburst tha t produces an optically thick wind 
dHumphrevs & Davidson! 1 19941) but it i s not well-established 
if L BVs do form pseudo -photospheres. iLeitherer et al.l ([l989) 
and Ide Koter et alJ d 19961) found based on NLTE modeling that 
the observed variability of the photospheric radius and effec- 
tive temperature cannot be due to the formation of a pseudo- 
photo sphere but must be induced by a sub-photospheric insta- 
bility. [Smith et alJ (|2004), on the other hand, argued that objects 
close to the bi-stability jump and close to the Eddington limit 
could form pseudo-photospheres. In our discussion we adopt the 
pseudo-photosphere hypothesis, i.e., that the apparent changes 
in R71's spectrum and photometry are caused in the outermost 
layers. However, we are aware that the discussion is still open on 
this issue. 

Estimations of R71's current effective temperature based 
on intrinsic colors are likely to fail becaus e LBVs and R71 
in particular show generally a UV excess (IWolf et al.l [l980t 
ILennon et al . 1993) and the reddening towards R71 is dis- 
put ed. For exampl e , the r elation between T e g and {B-V)o found 
by iNeugent et all d2010l) result s in J eff = 56 00-6000 K if 
we adopt (B-V) = 0.6 mag (http://www.aavso.org) and al- 
low for different extinction values. However, spectral features 
imply a higher effective temperature. To determine R71's ap- 
parent temperature we used the line-depth ratios described in 
iKovtvukh & Gorloval d2000t) . We determined a current effective 
temperature of r e fr = 6 700+400 K. This corresponds to an ap- 
parent spectral type o f F5 to F8 and a bolometric co rrection of 
BC = 0.08-0.13 mag dHumphrevs & McElrovll 19841) . 

We also compar ed the 2012 X-sh ooter spectra to Kurucz 
model atmospheres (ICastelli et al.lll997l) . We varied the metal- 
licity from [Fe/H] = to [Fe/H] = -0.5 (close to th e litera- 
ture value [Fe/H] lmc = -0.34, see iLuck et ail 1 19981) . the ef- 
fective surface gravity logoff = 0.0 to logoff = 2.0, and the 
extinction Ay - 0.15 to Ay = 0.63 mag. Two immediate prob- 
lems arise when comparing the X-shooter spectra to model at- 
mosp heres. First , LBV s have weaker Balmer jumps than normal 
stars. iGuo et ail (120051) suggested that this may be due to a de- 
crease of neutral hydrogen because of collisions in the dense 
wind, due to extra emission from the recombination of hydro- 
gen, or due to the more extended atmospheres in which the ve- 
locity field changes with the radius and the fi nes of the Lyman 
series could be spread out. IGuo & Lil (120071) showed that the 
continuum energy distribution around the Balmer jump is sen- 
sitive to the wind velocity law. Second, the observed spectral en- 
ergy distribution of LBVs in outburst may be a combination of 
two sources (the central star and the envelope) . IGuo et ail (12005) 
proposed that LBVs form a non-homogeneous envelope during 
outburst, which is not optically thick in all directions. The ob- 
served spectral energy distribution comes from the central star 
(mostly in the UV) as well as the optically thick part of the enve- 
lope (mostly in the optical wavelength regions). When compar- 
ing model atmospheres to the X-shooter spectra we thus focus 
on the continuum emission redwards of the Balmer jump. 
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Table 2. Parameters of R71. 



State 


Source 


my 
(mag) 


■D^LMC 

(mag) 


A v 
(mag) 


Teff 

(K) 


BC 
(mag) 


log (L/L Q ) b 


M Li 

(mag) 


Quiescence 


Wolfetal. (1981) 


10.88 


18.50 


0.15 


13 600 


-0.76 


5.30 


-8.53 




van Genderen (1982) 


10.89 


18.60 


0.15 


14000 


-0.95 


5.37 


-8.71 




van Genderen (1988) 


10.87 


18.60 


0.37 


14700 


-1.10 


5.53 


-9.10 




Lennon et al. (1993) 


10.88 


18.45 


0.63 


17 250 


-1.69 


5.86 


-9.94 


Outburst 1970-1977 


van Genderen (1982) 


9.79 


18.60 


0.15 


11000 


-0.38 


5.58 


-9.24 




van Genderen (1988) 


9.77 


18.60 


0.37 


8150 


-0.10 


5.57 


-9.20 


Outburst 2012 


this paper 


8.7'' 


18.50 rf 


0.15 f 


6 650 


0.11 


5.82 


-9.84 






8.7 C 


18.50 rf 


0.37 / 


6 650 


0.11 


5.91 


-10.06 






8.7 C 


18.50 rf 


0.63^ 


6 650 


0.11 


6.02 


-10.32 



References, (a) In the cases were only E(B-V) was given in the literature, we adopted the standard extinction l aw with i?y = Ay / E(B-V) = 
3.1 dHowarthlll983l): (b) v alues were calcul ated using £>Mlmc - 1 8.5 m ag; (c) http://ww w.aavso .org but see iGamen et al ] HoH); (d) e.g., 
lFreedmanetal.U2001l) : (e) I Wolf et al.ldl981h : (f) lvan G enderen et al. ( 1988); (g) lLennon et~aIUl993l) . 



Models with higher surface gravity \ogg have weaker 
Balmer and Paschen jumps. Because of the above described is- 
sues, w e make an apriori assumption about log g. iLennon et al.l 
(1993) found logg = 1.8 for R71 in quiescent state. The ef- 
fective surface gravity in 2012 is likely logg e ff,20t2 5= 0.5 based 
on the fact that R71's effective radius has increased by up to 
a factor of 5, see below, and that low effective gravities of 
logged = 0.5-1.0 were found for several LBVs in outburst. For 
example, Ster ken et al.l (1 19911) es timated log g e s - 0.55 for HD 
160529 and Sz eifert et al.1 (Il993l) found logger = .75 for R40. 
Both objects are similar to R71. iGuo et all (120051) determined 
log g e ff = 1 for R127 and logged = 0.5 for Rl 10. We thus adopt 
log geff ,2012 = 0.5 for R71, also on the grounds that for lower 
surface gravities the models become unstable against radiation 
pressure. 

Figure |4] shows de-reddened (Ay = 0.15,0.37,0.63 mag) 
2012 November X-shooter spectra and model atmospheres of 
different effective tempe ratures. The a n alytic formula for the 
mean extinctio n law by Cardelli etal.1 ( 11989b and R w = 3.1 
dHow arth 1983) were adopted to de-redden the X-shooter spec- 
tra. In each case the best fit leads to a too high Balmer dis- 
continuity, but see above. Unfortunately, we cannot discrimi- 
nate between the different values of Ay used in the literature. 
For higher extinction values the X-shooter spectra are well fitted 
by model atmospheres with higher effective temperatures. We 
find that model atmospheres with T e ff = 6 600+300 K fit best the 
spectral energy distribution of our X-shooter spectra, if we take 
into account the range of extinction values Ay = 0. 15-0.63 mag. 
The derived effective temperature from model atmospheres is 
very similar to the value from the line ratio measurements dis- 
cussed above. R7 1 's pseudo-photosphere has certainly reached 
much lower temperatures than during its 1970s outburst when 
the apparent temperature was not cooler t han T e ff ~ 8 150 K . 
However, we do not confirm the result by IGamen et al.l (1201 2b 
that R71 has currently an early G-supergiant spectrum. 

While there is uncertainty of R71's parameters during its qui- 
escent state because of uncertainties in the extinction Ay and the 
bolometric correction BC (see Table 0, its current temperature 
implies a small value for B C. We adopt a d istance modulus of 
DMhMC = 18.5 mag (e.g.. iFreedman et"al1l2001b . a visual ex- 
tinction of Ay = 0.15-0.63 mag, and a bolometric correction of 
BC = 0.11 mag, and find that R71's current visual brightness of 
wJv,20i2 = 8.7 mag results in Mb o i,20i2 ~ -9.8 mag to -10.3 mag. 
R71's luminosity of logL^ji t 2on/L = 5.8-6.0 is about its clas- 
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Fig. 4. Comparison of de-reddened X-shooter spectra from 2012 
November and Kurucz model atmospheres (with fixed param- 
eters [Fe/H] = -0.5, logg = 0.5). The continuum energy dis- 
tribution is well fitted for the extinction values Ay found in the 
literature. Higher extinction values require models with slightly 
higher effective temperatures. 



sical Eddington luminosity (log Lrji^em/Lq ~ 5.8-6.1). The pre- 
emption bolometric luminosity was A/b i,quiescence 8.5 mag to 

-9.9 mag if we adopt DMlmc = 18.5 mag instead of the varying 
distance moduli used in the literature^ The bolometric luminos- 
ity thus increased by 0.4-1.3 mag. 

The very low effective temperature of R71's pseudo- 
photosphere during the c urrent eruption b elow 7 500 K implies 
a very high wind density (lDavidsonlll987t) . LBVs in quiescence 
have mass loss rates comparable to normal supergiants of the 



4 The range of the pre-emption bolometric luminosity is larger be- 
cause of the greater uncertainty in the bolometric correction BC. 
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Fig. 5. Schematic upper HR-diagram (modified from 
iHumphrevs & Davidsonl The red solid curve is the 

upper luminosity boundary. Black solid lines represent LBV 
transitions. The positions of R71 during its quiescent state and 
its current eruption are indicated with blue solid curves, see 
text (R71's quiescent and 19 70s outburst locations placed by 
IHumphrevs & Davidsonl 1 19941 are also shown). The dashed 
blue curves show the transitions for the lowest and highest 
extinction values found in the literature; (a) Ay = 0.15 mag and 
(b) Ay = 0.63 mag. 



same temperature and luminosity (10 7 M Q yr 1 to 10 5 M Q 
), which increase 10-100 times during an outburst. R71's 



mass loss rate during its quiescent state was 3x10 7 M yr 1 

- 1 dWolfetal. 



and d uring its 1970s outburst 5xl0~ 5 M yr" 1 ([ Wolf et al. 
119811) . We use Figure 1 and Equation 4 in Davi dsonl d 19871) 
for a rough estimate of R71's current mass loss rate and find 
^R7i. 20i2 ~ 7xl0~ 5 M B yr -1 to 7x1 0~ 3 M yr -1 . Equation 

1 in IWolf et all (119801) results in a current mass loss rate of 
A^R7i,20i2 ~ 2xl0 -4 M yr . R71 thus increased its mass loss 
rate by a factor of about 1 000. The total mass of ejected ma- 
terial during the last 6 years may therefore be on the order of 
10~ 3 M . We find no evidence of an explosion. The spectra show 
no signatures of fast moving material and the terminal velocity is 
comparable to R71's quiescent state. Also, du ring its 1970s out - 
burst the terminal velocity was comparable (IWolf et al.lll98l1) . 
This indicates that the change in mass loss rate is caused by a 
tremendously increased wind density only. 

Th e effective radius increased by a facto r of about 5 from 81- 
95 R (IWolf et al.ll98lllLennon et al.ll993l) during the quiescent 
phase and is now on the order of 500 R Q . This radius estimation 
is based on the apparent temperature and luminosity during the 
eruption and may not be very meaningful because t he opacity 
in the pseudo-photosphere is mainly due to scattering (Davidson 
119871) . 

Figure |5]shows the location of R71 in quiescence and during 
its current eruptive state in the HR-diagram. The solid blue lines 
indicate the range of bolometric luminosities obtained for each 
state, if we adopt DMlmc = 18.5 mag and the different values 
for Ay and BC found in the literature, see Table [2] The dashed 
blue lines indicate the transitions for the two extreme cases. The 
transition labelled '(a)' is the transition for the lowest extinction 
value, Ay = 0.15 mag, and the transition labelled '(b)' is the 
transition for the highest extinction value Ay = 0.63 mag. R71 
could have potentially increased its bolometric luminosity by up 
to 1.3 mag. 

R71's visual magnitude variation on the order of Amy ~ 

2 mag implies that the star is experiencing a classical LBV 



outburst (Hum phreys & Davidsonlll 9941) . However, the effective 
temperature appears to be unusually low and the mass loss rate 
unusually high. Also, the common notion is that during classical 
LBV outbursts the bolometric luminosity stays constant, which 
is most likely not the case for R71's current eruption. With a 
growing number of observed LBV outbursts and with the real- 
ization that outbursts of the same object can differ considerably 
this commonly adopted classification scheme may have to be re- 
considered. 



5. Conclusion 

R71's current eruption differs from its 1970s classical LBV out- 
burst in several aspects and is challenging our view on LBVs. 
Both the visual lightcurve and the spectra show many differ- 
ences compared to its previous outburst. The visual light in- 
creased by Amy > 2 mag over the last seven years and the spec- 
trum indicates an unusually low temperature of R71's pseudo- 
photosphere, mimicking the photosphere of a late-F supergiant. 
Balmer and Fe II P Cyg profiles, normally observed during LBV 
outbursts and present during R71's previous 1970s outburst, are 
absent. Low-excitation forbidden emission lines and Fe I P Cyg- 
like profiles appear in late 2012. 

We found an apparent temperature of R71 in 2012 of 
r e ff,20i2 ~ 6 650 K. This is a decrease by more than 7 000 K com- 
pared to its quiescent state and more than 1 500 K cooler than 
during its 1970s outburst. The bolometric luminosity increased 

by 0.4-1.3 mag to Mb o i,20i2 9.8 mag to -10.3 mag and R71 

may have moved beyond the Humphreys-Davidson limit. R71's 
apparent radius increased by a factor of 5 to about 500 R Q . 

We estimated that R71 has a current mass loss rate on the 
order of Mki\,2Q\2 ~ 5xl0~ 4 M yr -1 , a factor of about 1 000 
higher than during its quiescent state and at least a factor of 10 
larger than during its 1970s outburst. The wind velocity of R71 
is on the order of 100-200 km s _I , similar to that of other LBVs 
while the pseudo-photosphere has a velocity on the order of a 
few km s" 1 . We find no spectral signatures of fast moving ma- 
terial, which may indicate an explosion. The large increase in 
mass loss rate during the eruption is therefore likely due to an 
increased wind density only, thus constraining instability mech- 
anisms. 

LBVs have been considered as extragalactic distance indi- 
cators because they belong to the most luminous stars (e.g. , 
iTammann & Sandagdfl%a ISandagdfl983l: fl984l) . IWolil (fl989h 
found a relation between the luminosity of LBVs and their out- 
burst amplitudes. However, with an increasing number of ob- 
served LBV outbursts it becomes clear that they can differ con- 
siderably - even outbursts of the same object. S Dor, e.g., dis- 
played an early F-type supergiant spectrum in 1999, which had 
never been reported before, while its visual luminosity only in- 
creased by 0.3 mag (Massev 2000). The current eruption of R71 
has a much larger amplitude than its previous outburst in the 
1970s. Thus LBVs appear to be unreliable as distance indica- 
tors. 

R71 is experiencing a much more powerful LBV eruption 
than during its 1970s outburst. We are in the unique position to 
witness such a rare event as it unfolds. Our X-shooter monitor- 
ing program on R7 1 will secure key data for our understanding 
of LBV outbursts and the properties of stars near the upper lumi- 
nosity limit. Observations in the infrared after the eruption has 
subsided are highly desirable because they will provide valuable 
information about the total amount of material expelled. 
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